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Abstract 
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' We consider the evolution of perturbed cosmological spacetime with multiple fluids and fields in Einstein 

^ , gravity. Equations are presented in gauge-ready forms, and are presented in various forms using the curvature 

^ ' ($ or <^x) and isocurvature (S'(ij) or 5(j>(ij)) perturbation variables in the general background with K and A. We 

\ clarify the conditions for conserved curvature and isocurvature perturbations in the large-scale limit. Evolutions 

of curvature perturbations in many different gauge conditions are analysed extensively. In the multi-field system 
we present a general solution to the linear order in slow-roll parameters. 



PACS numbers: 98.80.Cq, 98.80.Hw, 98.70.Vc, 04.62-hv 



> 

^ : 1 Introduction 
<^ . . . 

' Cosmological models with multiple components of underlying energy-momentum content are not only general but 
also necessarily appear in many cosmological situations including the early and the late evolution stages of our 
' world models. Assuming the gauge degrees of freedom are properly fixed, for the n component medium with fluids 
. or fields in Einstein gravity we anticipate a set of n coupled second-order differential equations for the scalar-type 
Oh' perturbation, a set of n (in general) coupled first-order differential equations for the vector-type perturbation, and 
one second-order differential equation for the tensor- type perturbation. The situations of the vector- and tensor- type 
5-H ' perturbations are rather trivial to handle, and, in fact, the situation is similar even in some classes of generalized 
c/3 . versions of gravity theories However, the scalar-type perturbation in such a multi-component system is 

' naturally more complicated, and often requires numerical methods. Not only the equations are complicated for the 
scalar-type perturbation, the behavior and the way of handling the equations depend on our choice of the gauge 
" condition which we need to make, fn either analytic or numerical handling of the scalar-type perturbations, it is 
. useful and is often necessary to examine the problems from the perspectives of different gauge conditions. 
' In the following, for the benefit of future studies, we would like to present various useful forms of the scalar- 

type perturbation equations considering arbitrary numbers of imperfect fluids and scalar fields with general mutual 
interactions among the components in Einstein gravity. Equations are presented in the gauge-ready forms which 
allow to take any of our favorite gauge conditions easily, and also allow easy translation of a solution derived in 
one gauge condition to the ones in other gauge conditions. We also present closed form equations in terms of the 
curvature and isocurvature perturbation variables in several useful forms, and derive general solutions if available. 
We recently have presented similar work on systems with multiple minimally coupled scalar fields . In the present 
work we present a complementary part now including the multiple fluids as well as the interactions among fluids 
and fields; we also extend by considering general K and A and the fluid formulation of the field perturbations 

The manuscript is mainly concerned with the formal aspects of perturbations in the multi-component situations. 
There are some advantages of investigating perturbations in a formal way keeping their most general forms. For 
example, we emphasize that the basic set of equations for a single component case remains valid even in the 
multi-component situations with additional equations for the components. Also, by including most general energy- 
momentum tensor in the fluids system we will show that the same equations remain valid for the mixture of fluids 
and fields, and sometimes even for the generalized gravity theories. By comparing the fluid system and field system 
we can clarify the similarity and difference between them. Although, we intentionally have kept the presentation 
formal, many new results are presented together with the known ones; see §o for a summary. Previous studies of 
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perturbations in the multi-component fluid system in a similar spirit can be found in § |, § 0, |. We start by 
introducing our notations, fundamental equations, and strategy in §|^. We set c = 1. 

2 Perturbed world model and basic equations 

We consider an arbitrary number of mutually interacting imperfect fluids in Einstein gravity. The gravitational 
field equation and the energy-momentum conservation equation are 

Gab = SnGTab ~ Agab, T^b = 0. (1) 
The energy-momentum tensor can be decomposed into the sum of the individual one as 

Tab ~ ^ r'(i)afc, (2) 



and the energy-momentum conservation gives 

I 

where (i) indicates the i-th component of n matters (including fluids and fields) with i, j,k, . . . = 1, 2, . . . , n. 

We consider the general perturbations in the FLRW world model. In the spatially homogeneous and isotropic 
background we are considering, to the linear order the three different types of perturbation decouple from each 
other and evolve independently. The presence of tensor-type anisotropic stress contributes as the source or sink of 
the graviational wave of the perturbed spacetime. The rotational perturbation of the individual fluid is described 
by the angular-momentum conservation where the vector-type anisotropic stress and the mutual interaction terms 
among fluids can work as the source or sink of the rotation of individual component. The most general situations 
of vector- and tensor-type perturbations are presented in ||. In the following we will consider the scalar-type 
perturbation only. 

A metric with the most general scalar-type perturbation in the FLRW background is 

ds^ = -a^ {l + 2a)dTj^ -2a^P,adr]dx°' gf^^ {I + 2ip) + 2j^^\p dx'^dx'^, (4) 

where a{r]) is the cosmic scale factor, a, (3, 7 and are spacetime dependent perturbed order variables. A 
vertical bar | indicates a covariant derivative based on g^^^ . Considering the FLRW background and the scalar-type 
perturbation, the energy-momentum tensor can be decomposed into the collective fluid quantities as 

T° = -(/2 + ^/i), r° = -i(M-Hp)^;,., T|^(j5+<5p)5^+(^lv"V;3 + i<5^)^(^), (5) 

where /i(= fi + J/i), p{= p + Sp), v, and n'-'^ are the energy density, the isotropic pressure, the velocity or flux 
related variable, and the anisotropic stress, respectively; k is the comoving wavenumber. In terms of the individual 
matter's fluid quantities we have 

I I 

and similarly for p, 6p, {fj, + p)v, and tt^'*-'; an overbar indicates the background order quantity, and will be ignored 
unless necessary. The interaction terms among fluids introduced in eq. (^) are decomposed as 

(5(i)0 = -a [(5(i)(l + a) + (5(5(,;)] , Q(i)a = J{i),a- (7) 

The equations for the background are: 

i^^ = ^M+^-4, ij=-47rG(^+p) + ^, (8) 
+ 2,H = Qi^), (9) 

where an overdot indicates the time derivative based on t defined as dt = adrj. K is the normalized spatial curvature, 
and H = d/a. The two equations in eq. (ph follow from Gq and — 3Gq components of the field equation in eq. 
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(|l|), respectively, and eq. (||) follows from eq. (^. By adding eq. (||) over components we have /i + 3H{fi + p) = 
which also follows from eq. (||). 

We present a complete set of equations describing the scalar-type perturbation without fixing the temporal 
gauge condition, i.e., in a gauge-ready form |8[ |], ||: 



K = 3 {-ip + Ha) + -^x, 

^ (fi + Hk — — 47rGdu, 

n X = 127rG-(^ -1- p)v, 

Q, rZ 

X + Hx- a~ ip^ SttG^tt^"), 

k + 2Hk + fsH - ^] a = 47rG {5^i + 3Sp) , 



5p(i) + 3H + 5p(,f 
1 



a imt) +P(j))w(j)J 



1 



a 



M(j) + Pii) 



^ (,) ^ 

- 3 p 



(10) 

(11) 

(12) 
(13) 
(14) 
(15) 

(16) 



where x = + aj). ip and k correspond to the shear, the three-space curvature and the perturbed expansion 
of the normal-frame vector field, respectively, see [Q, ^ for details. Vq and A are the covariant derivative and 
the Laplacian based on gf^. Equations (|l0[|lj) are: the definition of k, ADM energy constraint (Gg component 
of the field equation), momentum constraint (G° component), tracefree part of ADM propagation (G^ — ^S'^G'^ 
component), and Raychaudhuri equation (GIJ — Gq component), respectively. Equations ( [l5| , [l6| ) follow from the 
energy conservation (Tj-^-jq.^ — Q(i)o) and the momentum conservation (Tj.^^^.^ ~ Q{i)a), respectively. We note that 
eqs. ( p^Jl6 ) in the same forms are valid even in a class of generalized gravity theories considered in |2|; see eqs. 
(46,47) in 1^. Equations (p^p^ also remain valid in the generalized gravity where we can reinterprete the fluid 
quantities as the effective ones; see the Appendix B of ||^. 

By adding properly eqs. ( p"5|p^ ) over all components of the fluids, and using properties in eq. (^, we get 
equations for the collective fluid quantities as: 

/ h \ 

(17) 
(18) 

(19) 



6fi + 3H {6p, + 5p) = (/i + p) K — ?,Ha v 



[a^(/x+p)w]' _ k 

a'^'i^+p) a 



1 



p + p 



6p — — 



2P ~3K 



■ (s) 



3 fc2 



We introduce 

e(i) = Sp^,-) - c^,) = e = Sp- cl6fi, cl = p/ fi, 

W(^)=P(^)/p.(^), ■W = p/fl, 

where e(i) and e are the entropic perturbations. Equation ( [l5| ) can be written in another form as 



+ 3H (^C(^) - -I- 3iJ (1 + u;(,)) q(t)S(^) 



(-3i/e(,)+,SQ(,)), 



(20) 



where = iH{^Hi) + P(i))q(i), and we used W(^i) = -3H{cf.^ ~ u;(i))(l + w;(i))(l - Q(i)); %) = 

In the following we briefly summarize the gauge- ready strategy suggested by Bardeen B and elaborated in [Q ; 
see Q for a recent extension. Under the gauge transformation x°- — x°' + with ^* = (0 = rj) the perturbed 
metric and fluid quantities change as |^ : 

/ . ]f 

\ J a a 

Sfl = dn- 6il(^,) = - f^{^)t, Sp = Sp- p£,\ 6p(^i) = 6p(i) - _p(,)C*, 

(5Q(i) = (5(5(i) - = + (21) 
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As the temporal gauge fixing condition we can impose one condition in any of these temporally gauge dependent 
variables: a = (synchronous gauge), ip = (uniform-curvature gauge), k = (uniform-expansion gauge), x = 
(zero-shear gauge), v/k = (comoving gauge), Sfi = (uniform-density gauge), Sp = 0, W(i)/fc = 0, (5^(i) = 0, 
6p(^i) = 0, etc. By examining these we notice that, except for the synchronous gauge condition (which fixes a = 0), 
each of the gauge conditions fixes the temporal gauge mode completely. Thus, a variable in such a gauge condition 
uniquely corresponds to a gauge-invariant combination which combines the variable concerned and the variable 
used in the gauge condition. As examples, one can recognize the following combinations are gauge- invariant. 

6fiy EE 6fi~ ^fiv, ipy = ip-^v, ip^ = ip - Hx = -HXip, v^=v-^x, (22) 
etc. The original study of cosmological perturbation by Lifshitz M was made in the synchronous gauge. The 



zero-shear gauge and the comoving gauge were introduced by Harrison pQ l and Nariai 1 11 , respectively. The 
combination (p^ was first introduced by Lukash [ p2[ and Bardeen [ p^ . The gauge-invariant combinations ip^^ and 
are equivalent to p and v in the zero-shear gauge which takes x = as the gauge condition, and S/iy is the same as 
6fi in the comoving gauge, etc. In this way, we can systematically construct various gauge-invariant combinations 
for a given variable. Since there exist many gauge-invariant combinations even for a given variable, say for Sfj, we 
could have Sfiip, Sfiy, Sfiy^.^, etc., using our notation we can easily recognize which gauge- invariant combination we 
are using. Generally, we do not know the suitable gauge condition a priori. The proposal made in j^, [l], ^ is that 
we write the set of equation without fixing the (temporal) gauge condition and arrange the equations so that we 
can implement easily various fundamental gauge conditions: eqs. ( p^|l^ ) are arranged accordingly. We term it a 
gauge-ready approach. 

3 Fluid system 
3.1 Newtonian analog 

The following equations most closely resemble Newtonian hydrodynamic equations in the context of linear pertur- 
bation of the Friedmann world model: from eqs. (|ll|,^, eqs. (|l3,|l8|,^, eqs. (^Isj), eqs. (|io|,0|l|), and eq (^, 
respectively we have 

—p^ = AnGd^iy, (23) 



6(iy + SHS^y — 



k^-3K 



k 



2 



k 

-{fi+p)v-^ + 2Htt 



(24) 



k ( Spy 2k^ -3K n'-'^ \ , , 

2 

Hp^ = -'iTTG^ (/i + p)v^- SttGH^tt'^'^ , (26) 



fx 



a^ = -^^-87rG-7r(^). (27) 

Equations (p^psj) can be compared with the Poisson, the continuity, and the Euler (Navier-Stokes) equations in 
Newtonian theory, respectively. By setting p = = tt^'*-' and ignoring K = we recover the well known equations 
in the Newtonian perturbation theory with —p^, Spy and v-^ corresponding to the Newtonian potential, density 
and velocity perturbations, respectively. Notice the appearence of pressure terms in rather unexpected places in 
eqs. (^ 25) and none appearing in eq. (23). Although the zero-shear gauge and the comoving gauge were first 



studied by Harrison and Nariai, respectively ||10|, these equations using mixed gauge-invariant combinations 
were first presented by Bardeen in eqs. (4.3-4.8) of p3[. Newtonian correspondence of these fully general relativistic 
equations and variables was investigated in p4) . In Fig. |l| we present the behaviors of p-^^, (^(c)x) ^'^'^ '^v i^{c)v) 
in the case of three component system with the photon, baryon and the cold dark matter (c). We also compared 
the behavior with p, v (i'(c)) and d ((5(c)) in the uniform-expansion gauge. 

In we noted that in the sub-horizon scale in matter dominated era p, (5, and v in the uniform-expansion 
gauge all behave like p^, 6y, and which most closely resemble the corresponding Newtonian behavior, see also 
§84 of [h5|. We compared the behaviors in the uniform-expansion gauge with the variables with Newtonian behavior 
in Fig. lu 
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log(a/a„) 



log{a/a„) 



Figure 1: In the first figure, evolutions of (p^ (black, solid line), (red, solid line), (blue, dot and 

long-dash hne) are compared with f^. (black, dotted line), ^(^i^^^^ ^k (red, dotted hue), u„ (blue, dot and short-dash 
line). In the second figure, evolutions of ip^ (black, solid line), ^5(c)d (red, solid line), f(c)x (blue, dot and long-dash 
line) are compared with ip^ (black, dotted line), ■|5(c)k (red, dotted line), W(c)k (blue, dot and short-dash line). 
The vertical scale indicates the relative amplitude; we normalized = 1 in the early radiation era. We took 
K = ~ A, ilfco = 0.06, H — 65kni/secAIpc, and ignored the direct coupling between the baryon and the photon 
before recombination. Radiation-matter equality {teq) occurs around log(aeq/ao) ~ —4.2. We took the adiabatic 
initial condition in the radiation dominated era. As the scale we considered Xq = 2Trao/k = AirMpc. We indicate 
the horizon-crossing epochs of the given scale as ti where k/aH = 1, and t2 where X/Xr = 2naH/k ~ 1. 



From eqs. ( p^ , p7| , p^ we can derive 



1 + w 



- 



3K 



e l + w /a'^H'^ 7r(") 
- + 2- 



H 



k. 1^ + P 



2^ 
3 fi 



(28) 



The multi-component version of this equation will be presented later, see eq. (61). This equation is fully relativistic 
but most closely resembes the well known Newtonian density perturbation equation; only for pressureless fluid 
our comoving gauge coincides with the synchronous gauge. 6y is directly related to the density gradient variable 
introduced in the covariant approach] 1^; see §IV and the note added in [|l^. By using eq. (p^), the equation for 
ip^ can be derived simply from eq. (|2S), see eq. (|35|). 

In the multi-component situation, although eqs. ( ^ , ^5|) are still valid for the total fluid quantities, we can also 
derive the equations for individual component. From eqs. (|l^,^|2|) and eqs. (16jl^), respectively, we have 



3K 



2Ht: 



+127rG'-(/x+p)(^(j) +p(j)) {v^ - W(i)x) +<3(j)a«(i) 
Q(i) 



3HJ, 



Hv 

k 
a 



fJ-ii) +P{i) 



o^r^«' (a) 2P 



3K 



J, 



{i)x 



+ P(i) M(j) + P(i} 



(29) 



(30) 



By adding eqs. (^ 30) over the components properly we have eqs. ( |2^ , |25| ). Later we will notice that the above 
equations are not only valid for multiple fluids but also for the system with additional multiple fields. 



3.2 Curvature perturbations 

We introduce |l| 



which becomes Lpv, the Lukash variable, for vanishing K. Using eqs. (|2|j2^) we can show 



$ = 



4:'KG{p.+p)a \H 



2H^ a" 
II + p k'^ 



On the other hand, from eqs. 



3U3g) or eqs. (|22 



5^) we have 



H 



AiiGip + p) a? '^^ fi + p 



e TT 

3 



is) 



(31) 



(32) 



(33) 



Thus, in a pressureless case, if we could ignore e and tt*^^), $(x, t) — C(x) is an exact solution valid in general scale 
[p^ . Combining eqs. (^,^3|) we have closed form equations for $ and ip^: 



2„2 



(^+p)c 

IJ.+P 



r$ = - 



2„2 



{p+p)a^ 



e--vr 



2H^ 



if 



+ p) 



fx 



-4ttG ( e - -7r(^^ 



M + P 



II + p 



(34) 

(35) 



Equation ( |34| ) is valid for 0; for = we have eq. (|3|) instead. 

In the case of multiple ideal fluids without direct interactions among them, later we will see that the differences 
in the sound velocities cause nonvanishing e term which mixes the curvature perturbation $ with the isocurvature 
modes, see eq. (^6|). In §^we will see that a minimally coupled scalar field has a nonvanishing cntropic perturbation, 
see eq. ([76|). Thus, we can derive a similar equation for $ as in cq. (^) which is the form valid even in the multi- 
component case. 

Although introduced in a different way, a variable proportional to was first introduced by Field and Shepley 
in 1968 ||l^, see also [pOl. In the single component ideal fluid case, using w = z$ with z = +p/{csH) and 



d_ 



eq. (134|) can be written as 



v" + {cie - z"lz) i; = 0, 



(36) 



which is a well known equation in the cosmological perturbation, jl^, 19, [22| . This equation can be foimd in 
eq. (45) of |l^ by Field and Shepley. The variable H in eq. (43) of |lg| and in eq. (3.3) of |19 are proportional 
to V, and in eq. (122) of [0 is the same as — $0. 

The C, variable introduced in |j is the same as ips = ip + j^^^f^- From eqs. ( p^yi7| ) we can derive 



k He 
6a II + P 



(37) 



Using eqs. (^^ we have 

(2 - icDHifs 



•^s—'fS = 
0^ 



k-" 



He 
p+p 



M + P 



SttG 2k^-3K 7r(") 
7r(^) H 



M + P 



(38) 



We can also derive the closed form equation for ipg. From the definition of ips, evaluating it in the comoving gauge, 
and using eqs. (|2^,|3l|) we have 



(fisik^t) = $ + 



1 



a2 127rG(/i + p) 



(39) 



^ Ignorant of this rich history, and even forgetting our own contribution unfortunately, we have rediscovered this important variable 
* in pil. 
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where $ and Lp-^ are given above. For lPk, from the definition of and using eq. ( |11| ) and eq. we have 



1 



127rG'(Ai + p)a2 



1 



fc2 - 3X 



127rG'(/i + p)a2 



a77 



(40) 



Using the first expression in eq. (^8|) we can guess the equation for Lp^ . All the equations above are generally valid 
in the multi-component situation including fluids and fields. 

In the single component ideal fluid (later we will see that it also applies for multiple ideal fluids under an 
adiabatic condition, thus vanishing e) with ignorable c^fc^ term compared with z" / z term in eq. (36), thus valid in 
the super-sound-horizon scale, we have general solutions jl^, [Ti] 



$(fc,t) = C(/c)-d(fc)- 



C2i?2 



dt. 



<p^(k,t) = AnGCik)^ J' ''^^^/^ dt + d{k)^, 



(41) 
(42) 



where C(fc) and d{k) are the coefficients of the growing and decaying solutions, respectively. The coefficients are 
matched using eqs. ( ^2|j3^ ). In order to derive $ from ip^- and vice versa to proper order, we need to consider the 
next leading order terms in the growing solution of <I> and the decaying solution of ip-f^^; these follow trivially from 
eqs. (HH) as 



$ = C{k) i^l + k 



a{ii+p) 



dt dt 



a(M+-P) 



dt 



dt 



-d{k) 



A-kG J a^ifi +p) 
H f* a{p+p) 



+d(k)—ll + k^ 

a 



t 



dt, 
dt 



a(M+P) 



dt 



a^{fi+p) 



dt- 



dt 



a^{fi+p) 



dt 



(43) 



Equation (41) includes — limit; thus for a pressureless fluid we simply have $ = C, [jlSj. Solution for pg follows 
from eq. (|39|). Solution for ip,^ follows using the first expression in eq. (pO|). We emphasize that these solutions are 
valid in the super-sound-horizon (thus virtually in all scales in the matter dominated era). 

Notice that the decaying solutions of (ps, and ip^ (ignoring K in this case) are higher order in the large-scale 
expansion compared with the one of ip^. In the super- horizon scale we have 



<Pv ^ ips = ipK = G, 



(44) 



thus time independent, whereas ip^ evolves in time according to eq. (p2). As mentioned in §3.1, despite its complex 



behavior compared with the curvature variables in other gauges, ip-^ most closely resembles the behavior of the 
perturbed Newtonian gravitational potential . In the K = = A background dominated by an ideal fluid with 
constant w, in the super-sound horizon, the growing mode behaves as 



3 + 3u; 
5 + 3w;' 



C, 



(45) 



thus <£2: — f Ji 



and |(/9^ in the radiation and the matter dominated eras, respectively. We show the evolutions in 



Figs. Gl and y. The amplitudes of ip^, (p^ and p>^ drop as the scale comes inside the sound-horizon; compare Figs. 
^ and^. 

In the sub-horizon scale in matter dominated era we have the following behaviors: 

(i) ipv remains conserved. 

(ii) (ps approaches ^Sy. This is because the second term in eq. ( ^9| ) dominates the first, thus 'ps{= p>v + Jirfwj^v) — 

(in the sub-horizon K term is not important); it also follows from eq. (^). 

(iii) approaches (p^^. This follows from the second expression in eq. (^0|) and using eq. (p6|). 

We compared these behaviors in Fig. |^. Oscillatory behaviors of ip>v and others occur as the scale comes inside 
horizon during the radiation dominated era where the sound-horizon scale is comparable to the (visual) horizon 
scale. Since the above analytic solutions are valid for the super-sound-horizon we do not anticipate the above 
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Figure 2: The first figure shows evolutions of (^^ (black, solid hue), (p„ (red, dotted line), (blue, dot and 
short-dash line), ipv^c) (blue, dot and long-dash line). The same evolutions are reproduced in the second figure 
which also shows evolutions of (ps (red, solid line), fs^^^ (red, dotted line), ^d^^^^ (blue, dot and short-dash line), 
g'^(c)t)(c) (blue, dot and long-dash line). The conditions used are the same as in Fig. |l|. 



behaviors in (ii) and (iii) to be valid inside horizon in the radiation dominated era. The numerical results confirm 
( ii) and (iii) even in such a situation: Case (iii) is more understandable by examining the second expression in eq. 
( ^0|) and eq. (|2^) without using the large-scale solution. 

Although many of the results in this section were presented in our own previous work in |^ , |2^ , which were 
mainly concerned with the single component fluid case, we emphasize that by keeping the most general form of 
fluid quantities the equations are valid even in the context of the multiple fluids and fields as well. Figures || and ^ 
show how the conserved quantities in the single component case are affected by the multi-component nature of the 
more realistic situations. 



3.3 Isocurvature perturbations 



We emphasize that eqs. ( |23h40[ ) are generally valid for the multi-component situation. In the multi-component 
case e has additional contributions besides the entropic perturbation from the individual component. Following a 
pioneering study by Kodama and Sasaki, using eq. (|l^) we decompose Q 

e — Ci-cl + Cint, Cint = ^ e(fc), 
k 

ifJ-ik) + P(k)){fJ'{l) + P{1)) . 2 2 Nc , M(fe) +P(fc) 



_Y^/2 2Ne 1 (M(fc) +P(fe))(M(;) +^^(0)/ 2 2nc , M(fe) +P(fc) 2 r f^f.\ 

e,ei = 2^(C(^ - cJJa^W = 2 ^''(^^ " ''(O)^') + 2^ (46) 

where we introduced the following combinations 



^(.) - (47) 

S'(y) is gauge- invariant for (/(j) = 0, see eq. (^). Our definition differs from a gauge- invariant definition based on 
the comoving gauge {v/k = 0) originally introduced in coinciding only for qj^) = 0; for another gauge-invariant 
definition see ||25|| . 

Thus, we notice that the RHSs of eqs. (jsjjs^) vanish for the following situations: 

(i) single ideal fluid (e = = tt^*^), 

(ii) multiple pressureless ideal fluids, thus c^.^ — 0, 

(iii) multiple ideal fluids (e(i) = = 7r(i)) with the c^--, = cj-y implying c^--) — cj, thus Croi = 0. 
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Figure 3: The same as first figure of Fig. for scales Ao = AOnAIpc and AOOnAIpc, respectively. 



(iv) multiple ideal fluids with adiabatic perturbation, i.e., 5'(y) = for gjj) = 0; for gjj) 7^ it is not meaningful to 
set >S'(y) = without imposing the gauge condition. 

In these situations the large-scale asymptotic solutions in eqs. (41 4^) remain valid, thus $ is generally conserved 
in the super-sound-horizon scale considering general K, A and time-varying p(^). The case (iv) shows that the 
strictly adiabatic mode (S^if) = 0) is conserved in the large-scale limit for general number of fluid components. 

In terms of 5'(ij) and a gauge-invariant combination v^ij){= — from eqs. ( |20| , p^ ) we can derive Q 



(y) 



M (fc) + P(k) 
fl + p 



V(,j) - 3H [g(,)(l -I- cl^) - g(,)(l + cfjj) 
3 



{S(ik) 
fl + p 



Uk)) 



+ iH(q(^,) - - ^Hei^,^) + SQ(,j), (48) 



+ - 

a 



= 3iJ 

,1/2 2 NV^m+m^c- ,0 ^ , 2Ar^-3X (s) J 



(49) 



where e(ij), 7r|j'*j-), 5Q(^ij), and J(y) are defined similarly as S^ijy. 



°(i) 



l^{i)+P{i) m+PU) 



J, 



(0 



0-) 



A^W+m l^{j)+PU) 



(50) 



Although S'(ij), 6Q(^ij^ and J(y) are not gauge-invariant for nonvanishing (/(i), eqs. (^8|j4S|) are written in a gauge- 
ready form; these equations were derived in eqs. (35,36) of §. Under the comoving gauge condition based on 
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v/k = 0, removing Uv using eq. we can show that eqs. ( ^ , ^ ) reduce to (the corrected) eqs. (5.53,5.57) in 

0; typographical errors in Q] were corrected in eqs. (A. 37, A. 38) of |2^. We emphasize that, since eqs. ( |l5| , |l6| ) are 
valid in the context of generalized gravity theories considered in eqs. (|4^,^) in the same forms remain valid in 
the context of the generalized gravity theories as well. 

In the case of multiple field system, we have nonvanishing entropic perturbation of the individual component. 
Later we will show that the relative entropic perturbation e(ij) causes a curious situation of making eq. (E^) simple, 
and in fact making it identically satisfied, and eq. ( |48| ) alone gives the second order system; see eqs. ( [78[ - |80| ) . 

For (7(i) — 0, S[ij) is gauge- invariant. In such a case eqs. ( ^ , ^ 28 ]4^) form a closed set of equations for S(ij) 
and 5v] and v term in eq. (By) combine to give In this way, (5/i„ and Stij) may provide the curvature and 



(51) 



isocurvature perturbations, respectively. Using eqs. (|23|,|33[) we can express Sfi^ in terms $ and 9?^ as 



(5 Ail, 



fc2 - 3K 



1 - 3X/fc2 



+ e TT^''-' 

H 3 



where e is given in eq. ([46|), and eqs. ( p4y35| ) provide the equation for $ and ip^. In this way, $ (or ip-^) and S'(jj) 
may provide the curvature and isocurvature perturbations. Our use of these nomenclature is somewhat misleading 
in the sense that, actually, the conditions <S'(ij) = make the perturbation adiabatic, and the condition 4> = (or 
ip^ = 0) makes the perturbations isocurvature. [In the literature, often ip^ is chosen to indicate the curvature mode. 
(fi^ is simply related to Sfiy as in eqs. (|2^), and it is related to $ as in eq. (^.] We consider $ (or (p^^) and S(^ij) 
characterize the curvature (adiabatic) and isocurvature perturbations, respectively. For nonvanishing q(j_) we still 
have a freedom to impose a temporal gauge-condition depending on the situation, and we can use eqs. ( pi)[]l4| ) to 
replace the metric and matter variables in terms of the relevant curvature and isocurvature perturbation variables. 
For n-component fluids we can derive n-set of coupled second-order differential equations. 

In the case we could ignore the mutual interaction among the fluids for the background (thus, set (/(.j) = 0) eqs. 
(Mm give 



2- 2^4) + 4)) 



- -H{c- 



2 



-U) 



M(/c) +Pik) 



+ 



lfc2 



(c(i) + cl))S(^J) + (4) - 4)) E 



M(fc) +P{k) 



fi + p 



n + p 



(^S(^tk) 



5, 



Uk)) 



- ^72 ^ '^(.7)) 



(5 Ail 



fl + p 



2^ 



H {3He 



Ho) 



3 
2 

fc2 



)E 



fl + p 



3-ff(e(ifc) + e(jfe)) — (5(5 (jfc) — SQi^jk) 



2fc2 
3" 



Ji 



(y) 



fc2 "fe) 

where J/i^ is given in eq. (|5l|); equations for $ and ip^ are given in eqs. (| 



interaction among fluids (thus, SQi^i) 
( p2|) becomes 



which vanishes for 
(i) pressureless fluids, 



= J(i)) and ideal fluids (thus, e 



= TT 



a 



2(4) ^4)) 



Ai + p' 



(52) 

l^sl). Assuming no mutual 
the term in the RHS of eq. 



(53) 



(ii) fluids with c^^^ 



(iii) (5At„ = 0, 

(iv) most importantly, hy comparing eqs. (|5^j5l| , |28| ,^) or eqs. (^,^,|^,|6l) , we notice that, for K ^ Q, the RHS 
becomes negligible in the large-scale limit. Therefore, in the large-scale limit the curvature mode in the RHS of eq. 
( p2|) does not contribute to the isocurvature modes, thus, isocurvature modes decouple from the curvature one; this 
result was known by Kodama and Sasaki, see eq. (2.16) in |^; for another argument in a two-component situation, 
see Ip^ . This should be compared with the scalar fleld system where we expect less strong decoupling in general; 
see below eq. (§1) and H] . 
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In the cases of (i) and (ii), assuming ideal fluids (e(y) = = ti"!*])) and without interactions among fluids 
{^Q(ij) = = J(ij)), eq. ( p2|) becomes 



^(,,) + (2 - 3c2) HS, 



(ij) + C-s ^2 ^(ij) 



(54) 



where = c^^^. In such a case the RHS of eq. (|2^) vanishes, thus we have a general large-scale solution for 6^ 



= ^ E ^{k)&(k)v = - 3^)^ c(fc) y 



47rG 



(55) 



where the coefficients are matched using eqs. ( p3| , ^2| ). If w = W(^i-j = constant we have — w, in the large-scale 
limit the general solution of eq. becomes 



S(ij) = Cij{k) + Dij{k) 



dt 



t2-3-u) ' 



(56) 



We note that although equations in terms of 5y and S'(y) are decoupled in this case, the equations directly in terms 
of (5(i) are generally coupled even in the same situation, see eq. (|6^) as an example; similar situations in the scalar 
field system can be found below eq. (90). 

In the two-component ideal fluid system without mutual interactions, using S = <S'(i2), eq. (52) becomes 

2 _ C(2)(M(1) +P(1)) +C^i)(Ai(2) +P{2)) 



S + H{2- 'Sci) S 

„2 



(57) 



Various asymptotic solutions in the case of radiation (ideal fluid with p(^r) — kP(r)) plus dust (pressureless ideal 
fluid with = 0) system were studied thoroughly by Kodama and Sasaki in |5|, 

The adiabatic and the various isocurvature initial conditions in a realistic four component system with the cold 
dark matter, the massless neutrinos and the tighly coupled baryon and photon in the radiation dominated limit are 
presented in M. 



3.4 Component equations 

In this subsection we assume no interaction among fluids, thus gjj) = and <5(5(i) = = J(i) . We have 

LPS =(p- = (p+ — ^ -, 



where in the second step we used (/(^^ = 0. From eqs. ( pX)|jlq ) we can derive 



He 



M(j) + P{i) 



Using eqs. (Ohq) we can derive 



He 



+P{i) 



{2~3cl.)H'- 



3a2 







4)) 




e 


(i) 




^P{i) 



■is) 



2k'^-3K 
9 



^(0 



a Ai(i)+P(») 



(58) 



(59) 



(60) 



Using ips^.^ = eq. ( |60|) can be directly converted to a set of equations for (5(i)^, i.e., equations for (5(i) in 

the uniform-curvature gauge. We plot the evolution of ^pg^^^ for the cold dark matter (c) in Figs. |^ and |^. During 
the super-horizon scale fs^^.^ and cps are conserved like ipv and others. As they come inside the horizon they are 
dominated by the density parts; in the matter dominated era using eqs. ( ]60| , p8| ) we can show fs^^.^ and ips behave 

as ~ 3<^(c)ti — 3^(c)t)(^)- 
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It is also convenient to have equations for the individual component of density perturbation in a gauge-invariant 
combination with the velocity of the same component, i.e., in terms of ^(i)D(i) = ^(i) + 3iJ(a/fc)(l + By 
evaluating eqs. ( p^jTsHT^ ) in the = gauge condition we can derive 



H5. 



2 + 3^)-2u;(,)) 
-47rG(l + W{t))Y^ fHk) (1 + 3c^fe))(5(fe)„(^j 

k 

[('^W +Pw)e(fc) - im +Pik)) ew 

^(*) k 



5K e(,) 



3K 



f^(i) +Pii) I 3 ^(i) 



„2 it2_(s) 



(i) 



(61) 



with 



47rG(l + W(^i)) pi(fe) (1 + 3cff^))6(k)v^„ = 47rG(l + ^ Ai(fe) (1 + 3c^fc))5(fe)^,(^ 

k k 



M(fc)(l + 3c^fe))[(l + u;(fc)) 



(i) 3ifii;(i)^(i)„(.j 



(l + u)(,)) ((5(fe)^,(„ 



127rGiJ 

3Hw(^k)S(k)v^,,y 



E 



where we used 



which follows from eq. ( p^ using k^^.^ = k + (3ij — k"^ /a?){a/k)v(^i). We have 



(i) 



(62) 



(63) 



(64) 



Equations ( |6l| , |62| ) provide a closed set of second-order differential equations in terms of 's. This includes the 

case with one additional scalar field, see eq. ( [76| ) where — 6^,{i)s,jni^ — ~4>'fi)'^s<t>{i) ! fo^' multiple scalar fields 

we also have nonvanishing (5(5(i) in eq. ([77| ). In the single component situation Jfj-j^^.^ becomes 5^ which behaves 
similarly as in the Newtonian situation [ |14| . In the single component limit eq. ( |6l| ) reduces to eq. (p8|). Using eqs. 
( |lO[]l^ ) one can also derive a closed form equations in terms of 

For pressureless ideal fluids (wj^) = and e(j) = = 7r|^j^) eq. ( |6l| ) gives 



+ 2iJ -I- C?. 



47rG^/i(fc)(5(fc) = 0, 



(65) 



under the = gauge condition; thus the complete set of equations is made under mixed gauge conditions. Using 
eqs. (|To|-p^ we can show that eq. (^) is valid under the = or i> = gauge conditions as well. Thus we can 
regard eq. ( |65| ) is valid under any single comoving gauge condition of = or = 0; the tVc) = gauge is the 
same as the synchronous gauge without remaining gauge mode. In fact, one can derive eq. (pq) in the Newtonian 
context. 

In the case of ideal fluids with one component a pressureless fluid (w(c) = 0), eq. (^ij) for i — c becomes 



2^f<5(c)i,(,) - 47rG^/i(fe)(l + ^icl^,))S(k)v^, 



0. 



(66) 



Solutions in the case of radiation and dust were studied in 



4 Scalar field system 

Considering our previous study of the multiple scalar field system in Q the following can be considered as a sup- 
plement to p] . Besides considering general K and A the equations are presented in the context of fluid formulation 
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using fluid quantities; in fact, part of the fluid formulation was presented in Sec. 4.4.2 of By presenting the 
scalar field system together with the real fiuid system we can compare the similarity and difference of the two 
systems more easily. 

A system of minimally coupled scalar fields can be re-interpreted as a system of fiuids with special fluid quantities. 
For an arbitrary number of minimally coupled scalar fields we have 

Tab = ^ ((l){k),a<P(k)M - ^ffafc'/>(fc)'''(/>(fc),c j " Vgab, (67) 

where indicates i-th scalar field with i, j,k, . . . — 1, 2, . . . rt; V = V{(j)(^f^)) = V^(0(i), 4'{2}: • ■ • j 4>(n))- Additionally 
we have equations of motion for the scalar fields as 

^i^\-yi^ = ^^ (68) 

where = dV / {d(j)i^^)). 

We decompose the scalar fields as — (f)^i) + 5(l>{i) ■ To the background order, from eqs. (^,^,|8|) we have 

(fc) (fe) 
^(,) + iH^^r) + V(i) = 0. (70) 

Equations (|||,^,^ provide a complete set for the background evolution. In order to match with eq. (^) it is 
convenient to introduce fluid quantities of the individual field as 

+ P(i) = (t>li) , = -3i7 (Ai(i) + ) . (71) 

In this way we have = by definition. From eq. ( |7l| ) we have = (0(i) a-ndpjj) = (^(i) — . 
To the perturbed order, from eqs. (0J6^) we have 



k k 



k 



(;i + p)z; = -^0(,)J0(,), 7r(^)=0. (72) 



° k 



In the mixture of the fluids and fields, the effective fluid quantities of the scalar fields in eqs. (|69| , |7^ ) should be 
added to the total fluid quantities in eq. (||). Equation (68) gives 

5<P{t) + SH6<j)(^,) + —6(l){i) + 2Z y, {^{3)^^(0) = ^(i) ('« + ")+ (^20(,) + 3i7(/)(j) j a. (73) 
Equations ( p^[]T^|7^ , [73| ) provide a complete set. Under the gauge transformation we have 

Thus, in a single component case, the uniform-fleld gauge 5(j) = Q gives w = which is the comoving gauge condition 
In order to match with eqs. (|l5| , |l6| ) it is convenient to introduce fluid quantities of the individual fleld as 

'^^(i) = - + '^P(i) = </'(»)'^<^(i) = ""T^' '^(O = O-C^^) 



In this way, we can show e(i) in eq. (^^ becomes 

ew = (1 ('^6) 
With these definitions of fluid quantities eqs. ( |l5| , [l^ ) are satisfied with the following results 

^ E i^hk)^'t'{i) - (l>{i)^4>(k)^ , J{i) = 0, (77) 
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where we used eq. ( [73|) to derive eq. (15). From eq. (76) we can derive 



3H-V 



-(4) -4)) 



fi + p k 2 



1 



E 



/^(fc) +P(A;) 



(tk) 



(78) 



Thus, eq. ([i9| ) simplifies to become — 2i/w(y) = {k/a)S(ij'^ which is satisfied identically; these were found in 
[|l|. Introducing a gauge- invariant combination 



hi) kj) ' 



have = {k/a)S(f)(ij), and eq. ( ^8| ) leads to 



(y) 



A^(fc) +P(fc) 



{tk) + S4'{jk) 



-3H+ — 



= 3i/(4)-4))^, 

where, from eq. ( |77| ) we have 



(ifc) - ^(j)(fc) 



(79) 

(80) 
(81) 



Using c?.^ = 1 + 2V;(j)/(3i?(/)(j)) we have 



1 / 3 • • 



V(i)</'(j) 



(0 



(82) 



Using eq. (23), assuming K = 0, eq. ( p2[ ) becomes eq. (26) in iQ; in the present case it is valid for general K. 
From eqs. (19 



7^ we have e — {1 — c?g)5pL^ — 2J2k ^(fc)<5'/'(fc)ti- Using eqs. ( pi| , |74[ ) we have 



e = (1 - cl)Sfj,y 



H + p 



^y,(k)kk)4'li)^(t>(ki)- 



k,l 



Thus, eq. ( ^ gives 



3K 



1 - -iK/k"^ 



1 - 3(1 - C2)if/fc2 



M + P 
H 



M + P 



(fc) '?'(*:)'?>(/) 



(83) 



The RHS of eq. (|82| ) can be compared with the one of eq. (|52|), i.e., eq. (p3D. As mentioned below eq. (|53D, in the 
large-scale limit isocurvature modes decouple from the curvature one for the multiple fluid system; in terms of ip^ 
eq. ( |53| ) shows k'^ /a^ factor. However, it is different for the multiple field system; in terms of tp^ eq. (|8^) shows 
k^ jaf^ factor. Equation (HJ) remains valid with vanishing 7r'^*\ thus 



(/i + p)a^ 



(m + p)q^ 



/r2 



(p+p)a^ 



{fi+p)Hcl 



^V(k)kk)k(l)^'^(''l) 
k,l 



(84) 



where c\ = 1 — 3(1 — c'^)K/k'^. This equation can be compared with eq. (HJ). Since fc^ — > does not imply 
c^fc2 — > 0, contrary to the ideal fluid situation, we have to be careful in examining the large-scale limit for if 7^ 
case. From eq. ( p8| ) or from eq. (^TJ), using eq. (fz^), we can derive 



1+w 



a{fj, + p) \ H 



5v 



A-2 



fc2 - 3ii: 1 2 



fi fi + p 



i)09{ki)- 



(85) 
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Equations (|8^^5|) provide a complete set in terms of the curvature mode $ (or dv, or (p^ using eq. (p3|)) and 
isocurvature modes S4>(^ij) . Various analytic solutions for the curvature and the isocurvature modes in a system of 
scalar fields were studied in ||] assuming K = — A. Here we have considered general K and A in the background 
FLRW world model and presented the fluid formulation of the system. Literatures on the subject can be found in 
; some selected ones are in , and for recent additions, see ^ Q . 



4.1 Slow-roll: linear-order solutions 

We introduce the following slow-roll parameters (for different definitions, see Isl, |3^) 



H 



2,m 



(86) 



In §4 of we have presented various solutions of eqs. (82]8^) to the zero-th order slow-roll limit of e^, i.e., 
ignoring e^, thus ^(i) = — 3iJ0(i). In this case eqs. (82 84|) decouple from each other: see eqs. (29,30) in 0. 
We consider K = — K. The assisted inflation based on multiple fields each with an exponential potential |3q] 
has Ci — € — constant. Perturbations in the assisted inflation can be managed in simple way and were studied in 

In addition to zero-th order in e^, by further assuming zero-th order in e, we have 

^ (aVO' + =0, ^ + = 0, (87) 



where we used that V^(i)(fe) terms in eq. ( |82| ) are of the order e and e^. Thus, the curvature and the isocurvature 
modes are decoupled, and in the large-scale limit we have the general solutions 0] 



rdt, 64)(ij){k,t) = Ctj{k) - Dtj{k) 



ipy{k,t) =C{k)~ D{k) 

Notice that the non-transient solutions remain constant in time. From eqs. (|2^|7^) we have 

H 



rdt. 



l-L+p 



(88) 



(89) 



In the single component scalar field, equation for dcjy^p, which is Scf) in the uniform-curvature gauge, resembles most 
closely the scalar field equation in the given (i.e., without accompanying metric fluctustions) curved background 
[p7[. Equation (17) in |0] can be written using the slow- roll parameters as 



r . . 

)^ -I- 3HSij>(i)^ + -^54>(i)^ = ^ |^87rG(/)(j) (3 - e + e., + et) - '-^H^f^ik '^'/'(fc),^, 



(90) 



which is still exact. Thus, notice that although equations for and are decoupled to the zero-th order in 

the slow-roll parameters, equations for (50(i)<p are generally coupled to the same order. 

Now, to the linear order in the slow- roll parameters, and considering the non-transient solutions of eq. 
eqs. (11,11) give 



i(<.=*c)-4< 



2ff2 



(m + p)' 



k \ hi) hi) 



htlxAo) 



(91) 



(92) 



where superscripts (0) and (1) indicate the variables to the zero-th order and the linear order in the slow-roll 
parameters, respectively. Apparently, we can extend this perturbative procedure to the higher order in the slow-roll 
parameters. Therefore, to the linear order in the slow-roll parameters the general solutions in the large-scale limit 
become 



t 1 

a-' 
t 1 

—dt]s^}^dt 



\dt 

a-' 

-^dt 

a-' 



a^S'-^'^dt, 



a'S^'^dt. 



(93) 
(94) 
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We note that to the hnear order in the slow-roll parameter and considering the non-transient solution in the 
large-scale limit, the RHS of eq. (^2[), using eq. (^3|), vanishes. Thus, in such a case the isocurvature modes decouple 
from the curvature one, whereas the curvature mode is coupled with the isocurvature ones, see eqs. Studies 
in more general situation with nonlinear sigma type coupling are presented in keeping covariant forms in field 
space. 



5 Fluid and field system 

The formulations presented in ^ and ^ are valid for a mixture of arbitrary numbers of fluids and scalar fields. 
Thus, we will not rewrite the complete set of equations again. 

As an example, in this section we consider a two-component system of a fluid and a field with a general interaction 
term between them. As the Lagrangian we consider 



(95) 



is the fluid part Lagrangian which also depends on the scalar field as L(y) = L(^f^(Qmd, gab, (j^)- The variation 
with respect to gab leads to eq. (|l]) with an identification Tat = T(^f)ab + T(^<t,)ab- The variation with respect to 
leads to the equation of motion for cj) 

-V^ = -L(/),^ ^ r. (96) 

Thus, from eq. (||) we have Q(f)a — ~Q{4>)a and 

%)a;f, = Q(.f)a = L(^f)^^(t)^a (97) 

where we used i, j, . . . = 1, 2 = /, 0. As we decompose F = f + (5r, from eq. (0) we have 

Q(/) = f 0, 5Q^^f) = 5T^ + T{5^ - 0a), J(^-) = -T5^. (98) 

We have Q[^) = —Q(f) and similarly for SQ^^^ and 

By interpreting the fluid quantities properly, equations in §^ and ^ remain valid. For the background we set 
/X = -I- and similarly for p with 

m^l'^' + y^ p{t>)^\i^-y (99) 

Equation (^ remains valid. Equation (^ is valid for {i) — (/), and eq. ( |96| ) leads to 

+ 3i/</> + = -r. (100) 
For the perturbations we set 5^1 = Sfi(^f) + 5^[,f,) and similarly for (5p, {ji -f p)u, and t:^^\ with 



Equatio ns f|lO| - [l4|) are valid for collective fluid quantities. Equations (|l^,|l^) are valid for the fluid with (i) — (/), 
and eq. (|96|) leads to 

■■ • P • / ■■ -N 

S<j> + SHScI) + —S(t> + V4,4,S(t> ^ (f) {k + a) + [2(1) + 3H(f)j a - 5T. (102) 

After choosing the temporal gauge condition we can derive a fourth-order differential equation, or a set of 
coupled two second-order differential equations. As an example, we derive one form of such a set of equations based 
on a gauge condition = which is the comoving gauge based on the fluid velocity; in the following we sim 

ply write V(^f^ = as the gauge condition. Under our gauge condition V(^f^ = the perturbation variables are gauge 
invariant, i.e., 5(/)(= = and 5(j) = 54)^^^^, where, from eqs. ( p2| , [7^ , we have: 
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We assi 

in terms of (5(^) and 



an ideal fluid with w = P[f)/lJ'{f) ~ constant. Using eq. ( |16D for the fluid component we can express a. 



Using eq. (|l5| ) we can express k in terms of (5(^) and (50. Thus, eqs. (^,102) lead to 



^(j.) + (2 - 3u;) i?(5(^) + 



— ~6w(i? + i7^] -47rG^(/)(H-u;)(l + 3w) 



SttG 



2w0^5(^) + (1 + w) ( 20(50 - V^Sct) 






+ 3i7(50 + f ^ + ) (50 - [(1 - i«)0(5(/) -2w(^4> + 3H4>^ 6, 



3H + 167rG0^ - (50, 



(104) 



(1 + w^)M(/) 



ST 



l + w 
1 



l + w 



-ii.0-0fi:My+2 



3ff0 



r0^ 

r(50 " 
/^(/) 



(105) 



Although the RHS of eq. (104) contains (5(/) term, it can be replaced by using eq. (105). Notice that eqs. ( |104 , 105) 
are valid for general K and A. The LHS of eq. (104) is the familiar density perturbation in the comoving gauge 




10a) is also the familiar perturbed scalar field equation without the 



Thus, our w, 



(/) 



gauge choice allows simple equations for the 



derived in ||ll|, and the LHS of of eq^ 
metric perturbation; compare with eq. 
uncoupled parts of the system. 

This set of equations or eqs. (|6ll , |76| ) will be useful to handle the following situations: 
(i) Warm inflation scenario with nonvanishing interaction term, F and (5r, between the field and the radiation 
(u; = i)i|. 

(it) Time- varying cosmological constant simulated using the scalar field, often called a quintessence iPOl; in this case 



we may ignore the direct interaction between the fluid and the field F = = ST. An application of eqs. (104,105) 
to a system with an exponential type field potential is made in [[40[. 



6 Summary 

In this paper we have investigated aspects of scalar-type cosmological perturbation in the context of multiple 
numbers of mutually interacting imperfect fluids and minimally coupled scalar fields in Einstein gravity. Equations 
are presented using the curvature (<& or ip^) and isocurvature (5'(ij) or S(j)(^ij)) perturbation variables. It looks there 
exists no clear consensus about the curvature/isocurvature decomposition of perturbation in the literature. Since 
the equations are all coupled, such a decomposition may have clearer meaning in setting up the initial conditions. 
In and ^ we have shown that the equations can be meaningfully classified into the our definitions of the curvature 
($ or (/?^) and isocurvature (^(y) or (50(y)) modes; see below eq. (|5l|). We have shown that cither $ or can 
characterize the curvature mode depending on situations. 

We have presented the equations in a gauge-ready form where we have a freedom to choose one temporal gauge 
condition from the variables in eqs. ( pT| , [74[ ) or the linear combinations of them. In the general situation, eqs. 
( |lo[]l^ ) provide a complete set. In the case of scalar fields we also have eqs. (^,^2|). In terms of the curvature and 
isocurvature perturbations, eqs. (|3^j4^ , |49 ,^l|) with eqs. (0-|lJ) provide a complete set, and in the special case 
with multiple fields eqs. (|80|-p^ or eqs. (6l]7^ provide complete sets. In §|| we have shown how to consider the 
case with mutual interaction among fluids and fields. 

Some new results found in this paper are the following: 

(i) Gauge-ready formulation of the fluid-field system in Einstein gravity, §^ Hand |^. 

(ii) Several useful forms of the individual fluid equations, eqs. (^,^^ and § t3.4| . 

(iii) Analyses of various curvature perturbations in realistic situations of multi-component system, §3.2. 

(iv) Equations (^ , ^9| ) are valid even in a class of generalized gravity theories, see below eq. ( [l^ ) and eq. (50). 

(v) Extension of fluid formulation of the multiple field system, 

(vi) In the scalar field system the isocurvature modes are less decoupled from the curvature mode in the large-scale 
limit compared with the fluids system, see below eq. (p3|). The couplig term vanishes to the linear order in slow-roll 
expansion, see below eq. (|94|). 



(vii) Solutions for c^.^ = c^^-^ fluids, eqs. ( p4f|56| ). 

(viii) Solutions valid to the linear order in the slow-roll parameters, §4.1 
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The sets of equations in §H may be useful to analyze the evolution of structures in the world models with 
quintessence and in the warm inflation scenario. However, these sets of equations are the ones derived in certain 
gauge conditions, and when we encounter new problems we believe it is best to go back to the original set of equations 
in the gauge-ready form and see whether we could gather any new perspective from other gauge conditions as well: 
the set is in eqs. ( p^[]T^ , [73| ) . Corresponding set of equations in a gauge-ready form which is applicable in a wide class 
of generalized gravity theories including contributions from the kinetic components based on Boltzmann equations 
is presented in |^ . Specific applications will be made in future occasions. 
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